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Transparent media exhibiting anomalous dispersion have been of considerable interest since Wang, 
Kuzmich, and Dogariu [Nature 406, 277 (2000)] first observed light propagate with superluminal and 
negative group velocities without absorption. Here, we propose an atomic model exhibiting these 
properties, based on a generalization of amplification without inversion in a five-level DIGS system. 
The system consists of a A atom prepared as in standard electromagnetically induced transparency 
(EIT), with two additional metastable ground states coupled to the A atom ground states by two 
RF/microwave fields. We consider two configurations by which population is incoherently pumped 
into the ground states of the atom. Under appropriate circumstances, we predict a pair of new 
gain lines with tunable width, separation, and height. Between these lines, absorption vanishes but 
dispersion is large and anomalous. The system described here is a significant improvement over 
other proposals in the anomalous dispersion literature in that it permits additional coherent control 
over the spectral properties of the anomalous region, including a possible 10 4 -fold increase over the 
group delay observed by Wang, Kuzmich, and Dogariu. 



I. INTRODUCTION 

It has long been predicted that under certain condi- 
tions, optical media exhibit anomalous dispersion (ie. the 
index of refraction increases with increasing wavelength) 
[l|. In these cases, the group velocity of a pulse of light 
with appropriate frequency can be larger than the speed 
of light in a vacuum or negative @ ■ Negative group ve- 
locity implies that a smoothly varying pulse of light will 
appear to exit the dispersive medium before it enters [|| . 
This is possible because different spectral components 
of the pulse interfere strongly in the dispersive region 
and transform the leading edge of the pulse into an im- 
age of the pulse's peak. Thus it is only with regard to 
the overall profile of the pulse that superluminal propa- 
gation appears to occur; anomalous dispersion does not 
permit superluminal signalling and causality is not vio- 
lated [J II Ell- 
in 2000, Wang, Kuzmich, and Dogariu were first to 
observe light propagating with negative group velocities 
in a region of low absorption and minimal reshaping in 
a Rb vapor cell [1, Q . Their experiment was based on a 
system proposed by Steinberg and Chiao in which 
a A atom with the excited state coupled to one ground 
state by two far-detuned pump lasers is probed by a weak 
beam near the transition between the excited state and 
the second ground state. With appropriate initial con- 
ditions, this configuration leads to two narrowly spaced 
Raman gain lines in the probe's spectrum, with a re- 
gion of anomalous dispersion but low absorption between 
them. In 2003, Bigelow et al. reported superluminal (and 
subluminal) group velocities in a room temperature solid 
(specifically, an alexandrite crystal) [ll], Il3| ■ That same 



year, the claim that causality is maintained in cases of 
superluminal group velocities was supported experimen- 
tally by Stenner et al., who showed that the detection 
of a non-analytic point in an incident wave (representing 
new information) on the far side of a region of anoma- 
lous dispersion took at least as long as in the vacuum 
case hj. it is now well established that superluminal 
and negative group velocities are possible in otherwise 
transparent and non-interfering media without violating 
relativity or causality. 

In the decade since these first experiments, systems 
exhibiting anomalous dispersion have received consider- 
able attention in the literature, among theorists and ex- 
perimentalists. Some of these studies followed up on 
early predictions that strongly-driven two level atoms 
14 , [la . [lfl H3, EH and degenerate three level atoms 
13 . [201 ] can also exhibit regions of large anomalous dis- 



persion but low absorption (for a comparative study of 
systems of these sorts and additional references, see [2ljj ) . 
Others have involved coupling the excited states of a V 
atom (22j or the ground states of a A atom [23| and then 
using the coupling field to coherently control the sign of 
the dispersion, while others still have employed a mag- 
netic field to induce a Zeeman splitting, to similar effect 
[24 , |25| . Many of these subsequent studies have involved 
coherent modification of the Raman gain process pro- 
posed by Steinberg and Chiao by introducing a second 



31, 



321 or a second ex- 
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excited state [2J, |27J, |2J, |29j, |3j 
cited state and a third ground state 

In the current contribution, we present a novel pro- 
posal for producing a narrow, closely-spaced gain dou 
blet with an intermediate window of anomalous disper- 
sion, using dressed interacting ground states (DIGS) 
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DIGS systems [3o| are a generalization of the dou- 
ble dark resonances (sometimes called interacting dark 
resonances) introduced by Lukin et al. [36| , Double 
dark resonances occur when one ground state of a A 
atom is strongly coupled to the excited state (as in stan- 
dard EIT), and also coupled to a third ground state 
via an RF/microwave field. The result is a standard 
EIT absorption spectrum, with a new absorption peak 
located at zero probe detuning. New absorption nulls 
appear to cither side of this peak. These systems have 
been studied extensively and observed experimentally 
[H M, H, 11, SO, El, E3 • DIGS systems include a fourth 
ground state (to produce a five-level atom) coupled to the 
second A ground state by a second RF/microwave field. 
This additional coupling has been predicted to split the 
double dark resonance absorption peak into two sym- 
metric peaks located within the transparency window 
f34| . These peaks have widths and locations that are 
tunable by varying the RF/microwave Rabi frequencies, 
leading to additional control over the optical response of 
an atomic system. 

Here, we modify the model presented in [34| by in- 
troducing pumping terms. The pumping transforms the 
lines we describe in the earlier paper into the pair of 
gain lines. The system studied here has several benefits 
over previously considered examples exhibiting anoma- 
lous dispersion. For one, even in the cases were the 
Raman gain lines were shown to have tunable heights 
and signs, there was no independent control over their 
widths. Moreover, here the locations of the gain lines 
vary with the field strengths, not frequency or phase, 
and so they are easier to tune than a Raman system, 
where the location of the peaks depends on the laser fre- 
quencies. The lines we predict have tunable locations, 
widths, and heights (allowing the smooth change of the 
sign of both the absorption coefficient and dispersion), 
and so the current system permits a very narrow window 
with ultrahigh anomalous dispersion, or alternatively, a 
maximally broad window throughout which anomalous 
dispersion is present. We predict that in a cold Rb gas 
prepared in analogy to the experiment in Ref. [43| . one 
could observe a negative group velocity index two or- 
ders of magnitude larger than the largest yet observed 
of — 14000 [44J , and four orders of magnitude larger than 
the group velocity index reported by Wang, Kuzmich, 
and Dogariu [![. 

Since the gain described here is generated without ap- 
preciable population accumulating in the excited state, 
the phenomenon we predict can be thought of as a gener- 
alized example of amplification without inversion ( AWI) , 
which is often also referred to as las ing without inversion 
[Hill. Unlike standard AWI HI lH [11 , however , the 
current system requires only incoherent pumping to ex- 
hibit gain, as opposed to coherent pumping. We derive 
analytic results to describe the optical response as a func- 
tion of the incoherent pumping rate(s) in both open and 
closed pumping configurations. These analytic results are 
a central feature of the current paper, as previous presen- 




FIG. 1: Our five level model. An excited state \a) is coupled 
to two lower energy state doublets, {|6}, \b')} and {|c>, |c')}. 
fi M is the Rabi frequency of a strong control beam cou- 
pling I a) and |c); Qb and f2 c are the Rabi frequencies of 
two RF/microwave fields coupling the members of each of the 
doublets. We study the propagation of a weak probe beam 
with Rabi frequency f2 p <C f2 c , £lb, fi^ near resonance with the 
a) <-» 1 6) transition. The specifics of the decay and pumping 
schemes will be treated in sections MI Al and MI Bl 



tations of related systems (eg. double dark state systems 
with pumping) present exclusively numerical solutions. 

The remainder of the paper will be organized as fol- 
lows. In section [Hi we will review the base model that we 
consider in this paper. For further details on the model, 
see [34j. In section UTT1 we will solve the master equations 
for both open and closed pumping configurations. In 
section HVl we will derive the linear susceptibility for the 
probe beam. Here, we will present an explanation and 
analysis of the gain lines and the anomalous dispersion 
mentioned above. Section fVl will treat Doppler broaden- 
ing in the DIGS system. Finally, in section IVII we will 
offer some conclusions, including a discussion of possible 
experimental realizations. 



II. MODEL PRELIMINARIES 

Our base system is a five-level atom (see Fig. [1]) in 
which two sets of ground states, {\b), \b')} and {|c), |c')}, 
interact with a single excited state, \a). Transitions be- 
tween each of the ground state doublets are assumed 
to be dipole forbidden; the members of the hyper- 
fine ground state doublets, meanwhile, are coupled via 
RF/microwave fields (in what follows, we will refer to 
these as RF fields for simplicity) of frequencies ojb and 
lu c . These fields have Rabi frequencies SI5 and f2 c . One of 
the states in the {|c), |c')} manifold (without loss of gen- 
erality, |c)) is strongly coupled to \a) by a field with Rabi 
frequency f2 M . We study the optical response of a field 
near resonance with the \b) <-> \a) transition, with Rabi 
frequency fl p . We describe this system via the Hamilto- 
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for pai, from which we can derive the linear susceptibility 
for the probe beam in section [TV] 



H =- (w„|a)(a| + (uj b + u p )\b){b\ + {uy +u p - u b )\b'){b'\ 

+ (uo c + ^)|c)(c| + (u c > + - v c )\d){d\ - ft^\a)(c\ 
-fi 6 |&')<&|-fi c |e')(c| -np|o)(6|)+h.c. (1) 

The ~ here indicates that we have written the Hamilto- 
nian in a rotating basis defined to eliminate explicit time 
dependence. From the Hamiltonian we obtain the de- 
tunings for the control laser A M = oj a — uj c — v ^ and the 
RF fields, A;, = up — u>b — v\> and A c = w c i — uj c — v c . 
Most importantly, the detuning of the probe laser from 
the | b) <-* | a) transition 



(2) 



will be used to express the optical susceptibility of the 
probe field. 

For further details on this base system (including an 
explicit definition of the rotating frame) , see section II of 
[34j |. The model considered there functions as the start- 
ing point for the current paper; it is essentially identical 
to the model described here, except for the pumping and 
decay schemes, representations of which will be discussed 
presently. 

In what follows, we will be interested in the linear 
susceptibility expanded around the \a) <-* \b) transition, 
which can be written in terms of the corresponding den- 
sity matrix element, p a b- Since we are not yet con- 
sidering any nonunitary contributions, the density ma- 
trix can be described by the von Neumann equation, 
ihdp/dt = [H,p]. The equations of motion for the com- 
ponents of the density matrix follow immediately. They 
are given in Appendix A. In what follows, we will make 
reference to these in defining the models for each of the 
pumping configurations. 

In the bare basis, we will find that to first order in 
Qp/iln (with assumptions that we describe below) p a b de- 
pends on the solutions to six linked differential equations, 
which in general cannot be solved analytically. However, 
under the rotation defined in [3~i| , in which we diagonal- 
ize the {\b), \b')} and {|c), |c')} subspaces of the Hamil- 
tonian, these six linked equations decouple into two sys- 
tems of three equations each. The details of the change 
of basis are given in appendices B and C. We will distin- 
guished the diagonalizcd states by using capital letters, 
{\b),\b')} - {\B),\B')} and {\c),\d)} - {\C),\C')}. 
The transformation for each of the subspaces individu- 
ally is of the same form as the dressed states of a two 
level atom. 



III. SOLUTIONS TO THE MASTER EQUATION 

Here we present and solve the Linblad equation for 
the evolution of the density matrix in both open (section 
IIII Al) and closed (section [ill B[) pumping configurations. 
In each case, we are focused on finding a linear solution 



A. Open pumping configuration 

As mentioned in the introduction, the amplification 
presented here does not require population inversion. It 
does require, however, that population be distributed be- 
tween the ground states-specifically in {\b), \ b')} and|c'). 
The details of this requirement are most clearly mani- 
fest when we populate these states directly, via pumping 
from unspecified external states. This is the case we con- 
sider in the current section: the ground states \b) and \d) 
are pumped from external levels; likewise, decay occurs 
to external levels. Note that we would obtain essentially 
the same results if we added additional pumping to the 
levels \b') and |c) since fib distributes the population be- 
tween \b) and \b') while in the limit we are interested in 
(fl c -C 0^), any population in \c) is optically pumped 
into |c'). 

To model pumping from and decay to external states, 
we modify Eqs. I All by including terms to model direct 
incoherent pumping to \b) and |c'), at rates Tb and r c > 
respectively. 



: d fc 



l — rr— ~ ir-i 



at 3 

Relaxation terms are modeled by, 

■ d f>33 



j = b,c' 



at 
. apjk 
1 at 



~il3kP3k 



where jj is the decay from state \j) and jjk — + 

7fc) + Ijk i s the full off-diagonal relaxation term. 

Our strategy in what follows will be to assume that 
fl p <C fib, Qc, rb, TV «C ftp, and then use two applications 
of perturbation theory. First, we work at zeroth order in 
dp /flu, which essentially decouples the {\b), \b'}} mani- 
fold from |o), |c), and \d). The {\b), \b')} subspace can be 
solved exactly at this order. The {|a), |c), |c')} subspace 
is more complicated. However, under the approximations 
already described p c ' c i varies slowly relative to the other 
terms of the density matrix in this subspace, and so we 
can assume the other terms will follow it adiabatically. 
Then we can solve for the other terms to first order in 
Qc/Qft as a function of the steady state population of 
\d), Pc'c'i an d use these first order solutions to find a 
self-consistent solution for p^, c , valid to order (flc/fl^) 2 . 
Finally, we will move to the dressed basis introduced in 
section iHland defined in appendix B and solve for p a b to 
first order in Qp/Q^ using the zeroth order (in f2 p /f2 M ) 
solutions as source terms. 

We assume that the control and RF fields are on res- 
onance, A M = Ah = A c = 0. (A full general solution for 
the susceptibility with arbitrary nonzero detunings is in- 
cluded in appendix D. The derivation is identical to the 
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current case.) Then, to zeroth order in O p /0^, the equa- 
tions of motion for the {\b}, \ b'}} manifold can be written 



as 



■ dpbb . . _ fi 6 _ 

= tr b - iJbPbb + -^{Pbb' ~ Pb'b) (3a) 

-ilb'Pb'b' - -rrihb' - h'b) (3b) 



.dpb'i 



dt 
. dpbb' 
1 dt 



-Ibb'pbb' + ^-(pbb - Pb'b') (3c) 



ax 



These can be solved by writing them in the form .. 
—MX + A, which has a steady state solution of X 
M-^A. We find 



-st _ n{2 lb 'lbb' + ng) 
Pbb - 



2-lblb'lbb' + {lb + lb'W b 



Pb'b' — 



2iblb'7bb' + (76 + lb')^l 



Pll' = (pf'b)* 



-ir b jb'^b 



^Iblb'lbb' + (7b + lb')^l' 



(4a) 
(4b) 
(4c) 



To zeroth order in O p /f2^, our equations of motion 
for the {\a}, |c), |c')} subspace are, first for the diagonal 
terms, 



.£>Paa _ _. _ 

1> llapaa ^ \Pca Pa 



(5a) 



$ Pec - 11 / \ / \ 

l ~gf = ~ l lcPcc + —{Pea - Pac) + ^(Pcc' - Pc'c) 



.dp c > c ' . O c . . 

I Q t = !V - lie' Pc'c' ^\Pcc' - Pc'c) 

and for the off-diagonal terms, 
. dp ca 



(5b) 
(5c) 



dt 



. dpc'a 

■ dpe'e 



t^fc'cPc'c H~ 



-—(Pcc- 


. x n c . 

Paa) 2 Pc'a 


(5d) 


, *V 

+ —PCC 


2 Pea 


(5e) 


+ -^Pc'a ' 


+■ -ylPcc' - Pec) 


(5f) 



We assume that f2^,7 a 3> r2 c ,r c <,7 c ,7 c '. With these 
assumptions, it is clear that all of the equations of mo- 
tion, except the one governing p c ' c ', are dominated by 
the terms proportional to and 7 a . This is our justi- 
fication for the claim that p c > c > varies slowly relative to 
the other terms in this subspace, and that therefore p c / c / 
can be treated as a constant with respect to the other 
equations of motion. 

To zeroth order in fi c /fi M , we have two decoupled sys- 
tems of homogeneous equations. The steady state oc- 
curs only when p { ^ a = p c ? c = p a °] = pi^ = p c °J = 0. 
Physically, this makes sense, since the system is non- 
conservative and there is no external pumping to these 
levels/coherences. To first order in Q c /Q^, pi]} — Pec — 



Pea — again, since the equations of motion are un- 
changed at this order. The second system, however, now 
leads to nonzero steady state values. These are 



= 

Pc a Pc c 4 7c , c7c , a + Q2 



P ( V = 3*. , 

re'e He e 



^le'clc'a + 



(6) 
(7) 



This linear solution is sufficient to reproduce the effect 
that we are interested in. Taking the first order solutions 
as the steady states, we can solve self-consistently for 
pf, c ,. We find 



-St 

Pc'c' 



(8) 



We are now very nearly in a position to solve for p ab to 
first order in We move now to the dressed basis 

defined in appendix B. The assumptions already stated, 
A 6 = A c = 0, imply that ttf = Q 6 , ttf = Q c , and 
9b = C = 7r/4 (see appendix B for definitions of these 
terms). In order to handle decay analytically, we assume 
that jab — lab 1 - This is reasonable, supposing both ex- 
pressions will be dominated by j a > Ib,!^,!^^!^'- 
Moreover, we assume that ib ~ lb', l b ^ ~ l b i c , and 

lie' ~ 7b'c' so tnat we can take 7c6 ~ 7c6' = lc and 
-f c , b rts j c i b > — 7c/. In Appendix C we give full expres- 
sions for the decay and dephasing of the relevant den- 
sity matrix components in the dressed basis in terms of 
the i a b, lc, and 7c. These approximations may seem 
slightly arbitrary, but they permit both enough simpli- 
fication to solve the problem entirely, and yet contain 
enough nuance for an adequate analysis of the effects of 
decoherence on the phenomena we predict. Finally, we 
assume that ice' , lw , lc, lc *C Qp,, as would occur, say, 
in a cold atomic gas, or in a hot gas with a buffer gas 
present. 

We already have that p^ 1 } = 0; moreover, under these 
new assumptions, p^ a — — ^-p s J c i, which gives that (see 

PCa = -^rw^pf'c'- Meanwhile, 

nalization leaves these invariant, and so pbb ~ PB'B' ~ 

Taking these 



appendix B) p C 'a 

Pbb ~ pb'b' 



— and p BB , = {pb'b)* 



7b+7b' """far,- Kfnn, ■- ( Jb+Jb , )Q b ■ 

together, we find steady state solutions 
Pas - ^ (Vb {2ilc - 2A P - n b ) 



x (2z 7c , - 2A P - n b ) + n 2 c (p% - y B )) (9) 



and 



PaB' 



2Z- 



(<Pb (2«7c - 2A P + Q b ) 



x (2iic - 2A P + n 6 ) + n 2 c (pf, c , - y B )) (io) 
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where 



Z± = nl(2i-f C ' - 2A P t ^b) - (2i7afc - 2A P t fib) 



((*7c + »7C' - 2A P t fit) 2 - n; 



(11) 



and where we have defined tys = (for a gen- 

eral definition of ^3b, see appendix D). Meanwhile, to be 
consistent with the other approximations made thus far 
we should note that pf, c , simplifies to 



P% 



(12) 



27c'a^ + 7c ,fi2- 

p a b can be found simply from Eqs. [5] and [TO] via the 

2 



relation p ab = ^ (p a B - PaB>)- 



B. Closed pumping configuration 

We chose to present the open pumping configuration 
first because we feel it distills the important parts of the 
dynamics: as we will argue in section IIV1 population in 
two ground states, \b) and \d) is sufficient to produce am- 
plification of the probe beam. Thus the essential physics 
of the system is already present in Eqs. IMTT1 some read- 
ers may prefer to skip directly to section HVl now. How- 
ever, the theoretical literature has tended to focus on 
closed systems. Moreover, in the open pumping configu- 
ration we pump |c') directly, which begs the question of 
whether population will accumulate in \c') in the steady 
state if it is not directly pumped there. For completeness, 
we will now present a more theoretically natural case, in 
which atoms are pumped directly from |6) to \a), from 
which they decay to \b), |c), and |c ). In the appropriate 
limit (fib/A;, — > 0), the solution presented here is an an- 
alytic solution for the system described in, for instance, 

m 

In the closed pumping case, atoms are pumped from \b) 
to \a), and decay is internal to our 5 level subspace. We 
assume that the ground states are stable for the purposes 
of the current calculation. \a) is assumed to decay only 
to \b), \c), and |c'), with branching ratios a&, a c , and a c >, 
respectively (a& + a c + a c i — 1). Since a nonzero value 
of p^ c , is necessary, it is crucial that a c > ^ 0. The base 
equations of motion, Eqs. IAI1 now have contributions 



.dpaa 



-i{la + r )Paa + irp bb (13a) 



~ i{a-bla + r)p aa - irpbb (13b) 

at 



.dp ct 



lOi c -{aPaa 



.dpc'c' 

I a, ~ lOic'laPaa- 

at 



(13c) 
(13d) 



The definitions of the off-diagonal relaxation rates are 
unchanged from the open pumping case, except that now 
lb = 76' = 7c = 7c = 0. 



Our strategy here will be the same as in the open 
pumping case. To zeroth order in O p /f2„, however, the 
two subsystems of the previous case no longer decouple. 
But we can again make an observation about the time 
scales in the problem that will permit some simplifica- 
tion. As before, we assume that f2 M ^> f2b,f2 c S> r; 
moreover, we take j a to be sufficiently less that fi^ for 
it to be the case that aci a is about an order of mag- 
nitude smaller than f2 M . These considerations allow us 
to assume that pbb, pb'b'i Pbb'i and pc'c' vary slowly rel- 
ative to the other elements. We can thus assume that 
the rapidly varying ones follow these adiabatically. We 
again work perturbatively in f2 p /f2 M , and solve for p ac i 
and pc'c- Now, however, pbb and p c 'c' depend on p aa , 
and so to find a self-consistent second order solution for 
the slowly varying populations, we require a second order 
perturbative solution for the rapidly varying populations. 

We again work under the assumption that the control 
and RF field detunings vanish. First, note that it is clear 
from inspection of the relevant equations of motion in 
Eqs. lAll that in the steady state, p s b \ = pfi b ,- Meanwhile, 
we can solve for the rapidly varying terms perturbatively 
in flc/Cl^. To zeroth order, we again find two sets of 
decoupled equations. The ones describing p ac > and p cc ' 
are homogeneous and decoupled from the pumping at 
this order, and so they vanish. The second system now 
has an inhomogeneous term, p^, We find, 



dpec fiu 

l ~Qf = la c7aPaa + -^(Pca ~ Pac) 

dpea . - , ~ N 

l ~Qj~ = —llcaPca + — (Pcc - Paa) 



These have a steady state solution of 



(14b) 
(14c) 



Paa = Pbb 



r + (1 - a c )7a 



.(0) _ ~ / r ( 2a c7alca + Ica^l) 



Pbb 



P { ca = Pbb 



(r + (1 - a G ) 7a )7 ca fi^ 

ra c 7 a - ijca 
(r + (1 - a c )7 a )7 ca fi M/ 



(15a) 
(15b) 
(15c) 



To first order, the equations of motion for p aa , p cc , and 
Pea are unchanged, as the density matrix elements pro- 
portional to f2 c are zero to zeroth order. The equations 
for p c i a and p C ' C , meanwhile, become 



.dp C 'a 

'~dt~ 

.dpc'c 



a. 



a. 



-Hc'aPc'a + -TTPc'c 7T P 



(16a) 



2 2 

^Ic'cPc'c + ^Pc'a + %(# c , - ^) (16b) 
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These have a steady state solution of 



can write, 



.(i) = gc^cpjg - gy + 2i lc , c p^ a >) 

Pda ^Ic'clc'a + to* 

~(1) _ ^ C (2i7c'q(pic ) - g>cO + ^^P^) 

Pc ' c ~ 4 7c , c7c , + ^2 



(17a) 
(17b) 



In the open pumping case, it was only necessary to 
solve for the coherences to first order in f2 c /fi M in order 
to find the second order population p c i c i. In contrast, 
to find a fully self-consistent second order solution for 
p c 'c' and pbb in the closed pumping case, we require a 
second order solution for p aa and p cc . The equations of 
motion for p C ' a and p c 'c are unchanged at this order. The 
equations of motion for p aa , p CCl and p ca meanwhile are 
now 



dp 1 

S-7jp = ~i{la + r)p aa + irpf b - -Q^Pca - Pac) (18a) 



.dpcc 

.9 pea 



aPaa ~ n^c{p^l ~ P^J) + TT^Pca ~ Pac) 



(18b) 



-ijcaPca + \&n{Pcc ~ Paa) ~ Tp'cP^l ( 18 °) 



These are solved by 



~(2) = 2rfa + i(py c -p<£)n c 
Paa 2(r + (l-a c ha) 



(19a) 



p^ = 

rcc 



l 



2(r + (l-ach^ X ^ a ^ a + 2n > 11 

+ ztt c (2(r + 7-)7c«+nJ) (Pc'l-PcJ) 

+ n c n M (r + (1 - a c )7 a )Q5 QC < + ^„)) (19b) 

~(2) = -2iraclapf b + (r + la){p c ^ c ~ Pcc'R 

^ ca 2(r + (1 - a c ) 7a )f2 M 



47 ca 



(19c) 



The next step is to solve for and p c / c ; self- 
consistently, in terms of these steady state solutions. But 
first we can simplify the expressions already stated using 
our initial approximations that r ^ D, c ,^b ^ More- 
over, since the ground states are assumed not to decay, 
Ic'c = 7c'c- We assume that dephasing effects are small, 
and take 7 C < C <C r. Note that although we could in prin- 
ciple proceed through the next step of the calculation 
without making these assumptions, the expressions thus 
derived are unwieldy. Moreover, these assumptions will 
be necessary presently when we move to solve to first 
order in f2 p /f2^, and so it is expedient (and consistent, 
given our initial assumptions) to make them now. We 



~(2) 
Pda 



-in /2i7 c ' Q O c 

Pc'c = - ( ^2 I ■• 



_ (2) _ rQ^Pbb + 2'fc'a^lPc'c' 

Paa ~ (1 - acha^ 

X (rp bb (a c jalca + 0^) 



(20a) 
(20b) 

(20c) 



Pel 



(1 - achatll 

+ 2^f {2 lda nl + Jamaica " (1 " 0! c )fi£)) 



_ (2 s _ -~i(2~fc'aPc'c>Qc + racPbb^fJ 
Pca ~ (l-a c )Sll 



(20d) 
(20e) 



The equations of motion for pbb and p c ' c ' can now be 
written as 



.dpbb ., v_( 2 ) . ~ 

l-gf = H"b7a + r)p a > - irpbb 



(21a) 



^=^7.P-£? + ^-pS) (21b) 

Inserting the expressions for the second order rapidly 
varying terms and combining the two resulting equations, 
we find the condition that the steady state populations 
Pbb and pc'c' must satisfy. 



2j c ' a a b ( ' 
Pbb = Pc'd 



ra r ' V fi?, 



(22) 



This condition specifies a unique pair of populations 
when we impose the additional constraint that the sum 
of the populations must be 1. Then, 



Pc 



ra c nl 



Aabjc'a^c + rac^l 



(23) 



The other populations, meanwhile, can now be written 
as 



Paa ~ 



Pbb — Pb'b' ~ 

Pcc « 



2a 6 7 c / a O^ 



4a b j C 'a^c + ra c ^l 



(24a) 
(24b) 
(24c) 



Paa and p cc vanish because they are of order 

n 2 c r/(nl la ) < (n c /n M ) 3 « !• 

From here, the strategy is the same as in the open 
pumping case, and in fact, the solution in that case 
carries over wholesale. Eqs. [9l [TUl and QT] are gen- 
eral statements in terms of the steady state solutions 
for Pac, Pac, and the p a B, PaB' terms. In the semi- 
dressed basis introduced in section [Til we now have 
Pac = Pac = -^Tn^Pcc', just as in the open pump- 
ing case; meanwhile p B B = PB'B' = Pbb, since p b b = Pb'b' 



7 



and p w = 0. So if we take *p B = p B B — 4ab7e , a og+ra c n 2 
and consider the closed pumping solution for p c > c 'i then 
Eqs. [9lfTT1 hold, now as a function of the pumping from 
1 6) to \a), r. 



IV. LINEAR RESPONSE OF THE PUMPED 
DIGS SYSTEMS 



(1 + Re[x (1) (Ap)]) 1/2 . Considerations arising from the 
particular experimental set-up would determine density 
profile, cr(r), and number density of atoms, N, which 
are included to formally account for the particular op- 
tical thickness of the sample. D ab = e(a\e ■ x\b) is the 
dipole moment between |a) and \b), as a function of laser 
polarization, e. In our analysis we will focus on the di- 
mensionless reduced susceptibility 



The complex linear susceptibility, expanded about the 
| a) <-* \b) transition, is given by 



2a(r)ND ab . 

-p. Pab- 



(25) 



X^ determines both the absorption coefficient, a(A p ) = 
fcpIm[x^(A p )], and the index of refraction, n(A p ) w 



D 2 ab Na(v) 



A 1 ) = 



2~/ab . 

ir Pab - 



(26) 



where we have used the definition of the probe Rabi fre- 
quency, fl p — D ab £ p /K. This allows us to write the fol- 
lowing expression for the linear susceptibility as a func- 
tion of p c 'c' and ^Pb, which will vary depending on the 
pumping configuration. 



Jab 



(*Pu (2i lc - 2A P - fl b ) {2i lc > - 2A P - n 6 ) + n 2 c ($ c , - «p B )) 



ni(2i lc - - 2A P - n b ) - (2% lab - 2A P - n b ) f (i 7c + hc> - 2A P - n b f - siA 

(% B {2i lc - 2A P + Q b ) (2i lc , - 2A P + Q b ) + Q 2 C (pf, c , - <p B )) 
^2 (2i lc , - 2A p + il b ) - {2i lab - 2A p + n b ) ((i 7c + i lc , - 2A p + n b ) 2 - 



(27) 



As a check on this solution, note that in the limit that 
p c i c i — > (corresponding to r c / — > 0) and tyg — > 1/2 (cor- 
responding to normalized population beginning in \b) and 
| b')), we recover the solution presented in [341 ] . provided 
that 7c, 7c are small; this latter solution, meanwhile, 
reduces to the standard EIT solution in the limit that 
il b /A b , f2 c / A c — > 0. Note that it also amounts to an an- 
alytic solution to the double dark resonance system, [3(| , 
in the limit that fl b /A b — > 0. (See Appendix D for x^ 
for arbitrary A^, and A c ) 



In section IIV Al we will examine the imaginary part 
of Eq. showing how the system(s) solved above lead 
to gain lines in the appropriate limits; section [IV Bl will 
treat the real part of Eq. [27j including the anomalous 
dispersion present between the gain lines. This discus- 
sion requires us to estimate the values for the important 
variables in the problem. The Rabi frequencies of the 
coupling laser and RF fields are experimentally tunable 
over a large range. For the spontaneous emission rate, 
we take 7 a = 10 7 s _1 and will measure the Rabi frequen- 
cies and detunings in units of j ab . For the ground state 
dephasing rates the range of values are limited primar- 
ily by collision rates and therefore are temperature and 
density dependent. However for concreteness, we assume 
that jab, 7c, an d 7c are in the range 10 3 — 10 4 s _1 . 



A. Im(x): Gain lines 

In Figs. [2] and [3l we plot the imaginary part of the 
linear susceptibility, Eq. 1271 for various choices of pump- 
ing rates, in both open and closed cases. In both cases, 
we compare our analytic results with a direct numeri- 
cal solution of full density matrix equations of motion. 
We see that the additional levels manifest themselves as 
two tunable resonances located inside of the EIT trans- 
parency window. In general for arbitrary A;,, the new 
resonances are symmetrically located about A p = at 
the locations A p = ±ftf/2 = ± V / A 2 b + fif/2. For 
^^t,7a6 Ui,li c ,7c,7c, their shape is approximately 
Lorentzian, given by (for A M = A b = 0): 

( tfc/al + ichab \ 

y (A p t si b /2y + ( 7ab (ni/rtl + lc >ha b )) 2 ) 

(28) 

in the vicinity A p as ±fi&/2. Eq. [221 shows that we 
can expect absorption for Re(^Ps) > pf,^ and gain for 
pf, c , > Re(?Ps). Note that these conditions are both 
necessary and sufficient for absorption and gain, respec- 
tively, which implies first that no population need occupy 
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(a) Open Pumping (Analytic and Numeric) 




2 7ab 



(b) Open Pumping (Analytic and Numeric) 



0.2 




0.1 




-0.1 




-0.2 




-0.3 




-0.4 






(a) Closed Pumping (Analytic and Numeric) 
Imx 




2 7ab 



(b) Closed Pumping (Analytic and Numeric) 

Imx 



0.2 




0.1 




-0.1 




-0.2 




-0.3 






FIG. 2: (Color online) Here we compare analytic and nu- 
meric solutions for the imaginary part of the reduced suscep- 
tibility (corresponding to the absorption coefficient) in the 
open pumping case. The numeric solutions are represented 
by wide dashed lines, and the analytics by narrower solid 
lines. Agreement between them is quite good. In the top 
plot, we show the full spectrum in the cases where r c / = 
(absorption lines; blue) and r c > — .0077 a b (gain lines; green). 
In the bottom plot, we show a close up of one of the narrow 
features. The pumping parameters here, going from top to 
bottom, are r c i — (blue lines), r c / = .0027 a b (red lines), 
r c / = .0047 a b (brown lines), and r e / = .0077 a t (green lines). 
In all cases, r b = .00017^6, 0. b = D. c = .l-y ah , fi M = 27 a6 , 



and 76 = 7b' = 7 C 
dephasings vanish, 7: 



,7c 



10" 



ph 



"fab- 



We have assumed the 



0. 



the excited state in order for amplification to occur (thus, 
we find amplification without inversion) and second that 
coherent pumping is not necessary for this amplification 
to occur. Populations in the appropriate ground states 
alone are necessary. 

In the case of nonzero 7c, the widths of the features, 



labiSlll^l+lC'hab), 



(29) 



is the sum of the 'power broadening' term 7 f,f2j;/^u and 
the dephasing rate for \d) while the height is given by 



Ini[x (1) (±0 6 /2)] = 



rc'c' 



(30) 



The dependence on the population pf, c , is manifest in this 
expression: the height of the ultranarrow features varies 



FIG. 3: (Color online) These are analytic and numeric so- 
lutions for the imaginary part of the susceptibility in the 
closed pumping case. Again, the numeric solutions are rep- 
resented by wide dashed lines, and the analytics by narrower 
solid lines. The top plot shows the full spectrum as we vary 
r; the bottom plot show a close up of one of the narrow 
features. In both plots, moving from the top curve to the 
bottom, the parameters are r = (blue lines), r = .0057 a b 
(red lines), r = .0l7 a 6 (brown lines), and r = .047 a b (green 



lines). In all cases, fib 



fi c 



■ljab, fi p 



27 o6 , and 



ic = 7c 



7c 



76' b = 10 Jab- 



linearly with the difference in populations between \b) (or 
in general a function of the population of \b}) and \c'}. 
When the population in |c') becomes large, the sign of 
the Lorentzian reverses, and the absorption line becomes 
a gain line. 

The populations p c > c > and Re(*Ps) vary with the pump- 
ing rates for each of the two pumping configurations. In 
the open pumping case, the relationships are linear in the 
pumping rates and r c > ; the other states do not attain 
appreciable population. We find gain when 



> 



(2 7c , fi=+ 7c ,J)2) 



(31) 



The linear dependence on the pumping rate in the open 
pumping configuration permits additional control over 
the shapes of the lines. As can be seen in Eq. [301 the 
heights of the features are proportional to ilj?. In the 
case where f2 c becomes small (as is necessary to narrow 
the widths of the features), it is theoretically possible to 
counteract the corresponding suppression of the feature's 
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Closed Pumping Populations 




0.01 



0.02 



0.03 



0.04 7ab 



FIG. 4: (Color online) The analytic solutions (solid lines) 
for the populations in the closed pumping case, plotted with 
numerical solutions (dashed lines) to the full system of equa- 
tions. The lines with positive slope are the population of 
c') in each case; the other non-zero lines are the populations 
of |6), \b'). The dashed lines with vanishing population corre- 
spond to |a) and |c). In all cases, Qb — f2 c = ■ l'Jab, = ^ab, 
and 7c = 7c" = 7c C ' = Ib'b = 10" 4 7 ab . 



|a ), E =hw a 



\a_), E_=*a>—Jtf+a<l 



\B'), E B =heo i +m b /2 

\b), E b =tiw„ 

\b), E B =hco h -hn h i2 



FIG. 5: Energy level diagram that indicates transitions in- 
duced by the probe laser between the ground state manifold 
{\B), \B')} and the excited state manifold {|o+), \a~), |ao)}- 
Transitions to the gain state |oo) are indicated by dashed 
lines. The energy of the bare state \b) is also shown for refer- 
ence. 



height by increasing the pumping rates r& (for absorp- 
tion) or r c i (for gain). 

In the closed pumping configuration, meanwhile, the 
populations depend on a single parameter, r. Now we 
find gain for 



r > 



2a b j c < a fl 2 c 



(32) 



Here, it is useful to compare the analytic expressions for 
the populations to the populations found by direct nu- 
merical calculation (see Fig. 2]) . As we see, there is excel- 
lent agreement in the case of small r, as assumed. As r 
grows, both plots plateau, but there is a small deviation 
between the numerics and analytics. This arises because 
we assumed r <C Slb,f2 c . Indeed, it is surprising that 
agreement is acceptable for r > fif,, fi c in these plots. 

In the case where p c / c / becomes large, the presence 
of gain lines without population inversion can be un- 
derstood in terms of dressed states. Again assuming 
A c = and A M = 0, the eigenstates of the Hamilto- 
nian for the {\a), \c), \c')} subsystem, W = fruj a (\a)(a\ + 
\c)(c\ + \c'){c'\) - (h/2)(n iJl \a){c\ + n c \d){c\ + h.c), are 



l«+> 



l 

71 
i 

71 



(sin0|a) + |c) + cos0|c')) 
(sin0|a) - |c) + cos0|c')) 



cos ma) — sm o\c 



(33) 

(34) 
(35) 



where tan0 = f^/i7 c . The energies of the states \a±) 
are E± — hu a ± H^/f^ + ^1/2 while |a ) has energy 
Eq = hu a . Note that {\a), |c), \c')} is isomorphic to a A 



atom and so |arj) is a dark state of the sort familiar from 
STIRAP and coherent population trapping. To distin- 
guish it from the more familiar dark state composed of 
\b) and |c) responsible for EIT, we will call \ao) the gain 
state. Fig. [5] shows a schematic diagram of the energy 
levels of the dressed ground state manifold {\B), \B')}, 
which are coupled to all three states of the excited state 
manifold {|a + ), |o_), |oo)} y i & t ne probe. 

The eigenstate \a ) is decoupled from the control laser, 
so there will not be any destructive quantum interference 
in the probe absorption or emission for transitions to |ao). 
This explains why we find spectral lines at the locations 
corresponding to transitions from \B) and \B') to |do). 
Moreover, in the limit as il c <§; f2 M , we find sin^ ^> cos6>, 
and so |ao) ~ \c'). Thus in the dressed state basis, the 
population in \c') corresponds to large population in the 
gain state. The gain state energy remains huj a , how- 
ever, even in this limit, which means that although there 
is no population inversion in the bare state basis, there 
is large inversion in this dressed state basis. Whereas 
when the population in \c') is small, transitions from the 
{\B), \B')} manifold to |arj) lead to absorption resonances 
at uj a — u>b,b'- These change to gain lines as the popula- 
tion in \d) ~ | do) increases. 

It is worth emphasizing that the gain lines described 
here arise from different physical process than, say, stan- 
dard driven two level pump-probe spectroscopy. There, 
gain arises from the exchange of light quanta between the 
strong pump field and the weak probe. The dressed states 
used to describe and explain pump-probe spectroscopy 
necessarily involve the probed transition. Here, the ba- 
sis in which the gain state emerges does not include the 
probed transition. In the present system, the gain is due 
to stimulated emission from a metastable quantum state, 
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|a ), to a decoupled (to zeroth order in fi p /fi M ) ground 
state. It is the coherent preparation of the {\a), \ c), |c')} 
manifold that permits emission from this state (essen- 
tially |c'), in the limits we have considered) even though 
the population of the bare excited state is negligible. 

The mechanism here also differs from textbook ampli- 
fication without inversion, as presented for instance by 
Scully and Zubairy [47|]. There gain occurs because the 
coherence between the ground states of the A atom (\b) 
and |c), here) cancels absorption, leaving only stimulated 
emission. Thus standard inversionless amplification sys- 
tem requires coherent pumping. This permits gain even 
in the presence of small population of the excited state, 
| a). But crucially, some population in the excited state is 
necessary. In the present system, \a) may have zero pop- 
ulation, so long as |c') is populated. We note that there 
are other approaches to amplification without inversion 
(which is ultimately a term that describes systems ex- 
hibiting a certain property, rather than the name of a 
particular class of systems) that are closer in spirit to 
the current system 46] . But it is nonetheless worthwhile 
to distinguish the DIGS approach from the "standard" 
amplification without inversion, since their mutual con- 
nection to A type atoms may produce confusion. 
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B. Re(x): Anomalous dispersion 

When the pumping rate is such that the gain peaks are 
large we find a region of anomalous dispersion but low 
absorption. This result should be unsurprising, as the 
Kramers-Kronig relations guarantee that any pair of suf- 
ficiently narrow and tall gain lines will give rise to anoma- 
lous dispersion, given the analyticity of the susceptibility 
[Tfj| . The dispersion is proportional to the first deriva- 
tive of the linear susceptibility. In the present case, the 
heights, widths, and separation of the peaks, and thus 
the magnitude and range of the anomalous dispersion, 
are all tunable by varying r (or rt> and r c '), fib, and Q c . 
The real part of the susceptibility is shown for several 
values of the pumping rate in Fig. [5] Our definition of 
the probe detuning, 

Ap = UJ a — LOb — Vp, 

implies that the dispersion is anomalous when the slope 
of the real part of \ is positive. 

Near the bare transition frequency, A p = 0, the real 
part of the susceptibility is approximately linear. To sec- 
ond order in Qb/fl^ and f2 c /f2 M (the order to which the 
solution for the susceptibility is valid), we find, 



Mx (L> ) 



A p . (36) 



See Fig. [7] for a comparison of this linear solution and the 
full real part of the susceptibility for several populations. 
This relation can be rewritten to reflect the relative pop- 
ulations necessary for dispersion to become anomalous. 



FIG. 6: (Color online) Here we compare our analytic and 
numeric solutions for various pumping rates in the open con- 
figuration (a) and closed configuration (b) in the vicinity of 
the new features. In both cases, we see that increasing the 
pumping (or, increasing pumping to |c'}) results in a change 
from normal dispersion (negative slope, on our sign conven- 
tions) to anomalous dispersion (positive slope) in the region 
around zero detuning. Once again, the analytic solutions are 
represented by narrow solid lines, and the numerics by wide 
dashed lines. The parameters in both cases are the same as 
in the close-ups in Figs. r2](b) and0(b). In particular, for (a) 
we have (from most negative slope to most positive), r c i = 
(blue lines), r c i = .0027 a f, (red lines), r c i = .0047 a f, (brown 
lines), and r c i = .0077 a 6 (green lines). For (b), we have r — 
(blue lines), r = .0027 a i, (red lines), r = .0047 a (, (brown lines), 
and r — .0077 a {, (green lines). 



The dispersion is anomalous just in case 

— > 1 + — b -. 

Re(<p B ) T ni 



(37) 



Eq. [37] can be rewritten as a constraint on Tb and r c / in 
the opening pumping case and r in the closed pumping 
case by substituting the expressions for the populations 
derived in sections IIII Al and IIII Bl In the open pumping 
case, the constraint on r c > jrb is given by 



r_d_ 

n 



> 



(27^ + 7^)^2 + ^2) 



(38) 



whereas the corresponding constraint on r in the closed 
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FIG. 7: (Color online) Here we compare our linear approx- 
imation of the real part of the suseptibility, Eq. 1361 (solid 
lines) with the general analytic solution, Eq. [27] (dashed 
lines) for various populations. For complete generality, and 
because the purpose is to show agreement between the lin- 
ear approximation and the full susceptibility, we consider the 
populations directly, here, rather than limiting ourselves to 
one pumping scheme or another. In order from most negative 
to most positive slopes, we have = 1 and p c / c i = (blue 
lines); = -25 and p c r c / — .5 (red lines); and ?P_b = -1 and 
Pdc! = .8 (brown lines). In all plots, yc = 7c = lab x 1CP 4 , 

= ^Jab, &b = Q c — Jab/W. 



FIG. 8: (Color online) Here we compare our linear approx- 
imation of the real part of the suseptibility, Eq. 1361 with 
the general analytic solution, Eq. 1271 as dephasing increases. 
The (blue) dotted line has yc = 7c = lab X 10~ 4 ; the (red) 
dashed line has yc = 7c — Jab x 10 -3 ; the (yellow) broad- 
dashed line has yc = 7c — 7ab/100. The solid (green) line 
is the linear approximation. Note that although the suscep- 
tibility is depressed in the vicinity of the new resonances as 
the dephasing increases, the dispersion in the linear regime 
does not change. In all plots, we take the generalized popula- 
tions p c i c i — .8 and Re(^Ps) = .1, while other parameters are 
fi M = 27a&, f^b = fi c = 7 a b/10. 



pumping case is given by, 



r > 



2a b j c > a (Ql+n 2 c ) 



fl 2 



(39) 



It is tempting to consider the limit that fi^/f^ 2 . — ► 0. 
However, Eqs. [M] and [37] are not valid in this limit. 
For one, taking — ► makes sense in the context of 
Eq. [37] but not Eq. [36] More importantly, the as- 
sumption of linearity only holds when there is a trans- 
parency window between the two peaks, or roughly when 
f2g/fi? <C SlbHab, which corresponds to the case where 
the separation of the peaks is larger than their widths. 
It follows that, while inversion of the populations of 
\<k>) ~ l c ') an d \b) is sufficient to produce gain lines, the 
constraint on the population of |arj) necessary for anoma- 
lous dispersion is more stringent. This result can be seen 
clearly in comparing Eqs. l38l and [3"51 to the corresponding 
expressions in section HV A[ Eqs. |3"T1 and l3"2l The pump- 
ing rate necessary for producing anomalous dispersion 
increases as the square of the width of the window be- 
tween the peaks, supposing that the widths of the peaks 
are held constant. Thus, though it is in principle possible 
to produce wide spectral regions of anomalous dispersion, 
there is a practical barrier imposed by how rapidly one 
can pump atoms into |c'). 

A distinctive and important feature of Eq. [32] is that 
the dispersion does not depend on the dcphasings yc 
and yc ■ This result is not a relic of the approximations 
that went into deriving the linearized equation, as can 
be seen in Fig. [5] The effect of dephasing between |c') 



and {\b), \b'}} is to destroy the coherences responsible for 
producing the new narrow resonances. But even as the 
resonances shrink, the region between the peaks near zero 
detuning remains unchanged, until yc> 3> ^ c and the 
resonances vanish altogether. 

The values of the anomalous dispersion that we predict 
have dramatic consequences for group velocity, which in 
the current context can be written v g = c/n g , where n g 
is the group velocity index, 

ng ~ n 2n dA p ■ 

n = (1 + Re(x^ 1 ))) 1 / 2 is the index of refraction, as de- 
fined above. The absolute value of the dispersion in the 
anomalous regime, for large population |c'), is compara- 
ble to the magnitude of EIT dispersion, as can be seen by 
examining Eq. [36] In the case where *Ps = 1/2, the dis- 
persion becomes approximately that of "standard EIT," 
which we take to be EIT in an identical system, without 
the two additional fields/levels we have introduced here. 



9Re(x«) 4 7o6 He(5p B )(n2 + n2) 



dAr. 



1 w 



(40) 



In the opposite limit, of very large pumping, p c i c i — > 1 
and we find 



1 v> 



(41) 



For fl c = fif,, Eq. |4"T1 reduces to the expression for EIT 
dispersion, with opposite sign. This means that if one 
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FIG. 9: The ratio r° IGS /r| IT as a function of the closed 
pumping rate, rj^ ab . We see that for r > 7 a b/10, the time 
delay becomes negative; for r > 37^/10, the magnitude of the 
negative time delay comes within a factor of 2 of the EIT time 
delay, permitting much faster light than previously observed. 
Here Q b = Q c = 0.l7 a b 



begins with a suitable EIT system and then introduces 
the additional couplings and pump processes, one can 
generate negative group velocities of the same magnitude 
as the ultraslow light observed by [H(J, and |5l| . 

Negative group velocities can best be understood in 
terms of the group delay, 

r d = e(i/v g - 1/C) 

where I is the sample thickness. To see the point made 
at the end of the last paragraph most clearly, suppose 
that a given EIT system exhibits a group delay of rf lT . 
Then, for index of refraction n rs 1 (as we have here) we 
can expect a group delay in the same system prepared 
with the additional DIGS couplings of, 



.digs = EIT Bc(VB)(nl + n 2 c )-n 2 c p. 

d ~d 2 



(42) 



We plot t^ 1gs / t^ igs as a function of the closed pumping 
rate in Fig. [9l 

To be perfectly concrete, take as a sample system the 
87 Rb vapor cell prepared in (43|, modified to include the 
additional couplings (see section IVIjl . We can assume 
that 0(, f» f2 c , so that the details of the strengths of the 
couplings are irrelevant. They observe group velocities of 
90m/s through their 2.5cm long sample, corresponding to 
a group delay of .26ms. In the corresponding DIGS sys- 
tem, with pumping such that p c i c i rs .8 and tys ~ -1 (this 
corresponds to a closed pumping rate of approximately 
r = .047 ab for Cl b « tt c « .l7 ab and S1 M = 2^ ab ), we 
should expect a group velocity of -150m/s, correspond- 
ing to a group delay of -.156ms and a group velocity 
index of n g = — 2 x 10 6 . These numbers represent an 
improvement of four orders of magnitude over Ref. Q, 
who found a group velocity index o f n„ ss —310, and of 
two orders of magnitude over Ref. [44j who observed a 



por system. This latter result is the largest superluminal 
group velocity yet observed directly. 



V. DOPPLER BROADENING 

Doppler broadening is an important experimental con- 
straint on the current system. Thus far, we have dis- 
regarded Doppler broadening in our treatments of DIGS 
systems, and so we will focus on it here. For a single pho- 
ton with wave vector k incident on an atom moving with 
velocity v, the Doppler effect shifts the atomic transition 
frequency loq by 

LO]J = UJq + k ■ V. 

The "Doppler width," then, is given by the width of the 
atomic velocity distribution od = ({^d-^o) 2 ) = uo<j v /c. 
For a Maxwell-Boltzmann distribution a v is given by 
(FWHM) 



<r tl =2W21n2— 
V m 

where fc^ is the Boltzmann constant, T is the tempera- 
ture, and m is the mass of the atoms. 

From these considerations, we conclude that Doppler 
shifts will be unimportant for the RF fields. Working, for 
instance, at room temperature (T = 300K), in a gas of 
Rubidium atoms, we would find that <r v « 400m/s. For 
an RF field, we can take luq ~ 100MHz, which gives an 
estimation of the Doppler width at ctd ~ 100Hz. Thus 
for the RF fields, the effect of Doppler broadening will 
be much less than the homogeneous broadening from de- 
phasing (which we have thus far estimated in our plots 
at 7 b , 7b/, 7c, 7c ~ 10~ 4 7 a fc ~ 10 3 Hz). In the same gas, 
however, taking ujq ~ lPHz — a characteristic frequency 
of light — we find that the Doppler width will be approx- 
imately an ~ 1GHz3> "fab- Thus Doppler broadening 
will dominate the optical transitions {\b), \b')} <-> \a) 
and {|c), |c')} \a). In a laser cooled system, where 
T w l^iK— lOO^iK, the Doppler width would be reduced 
to approximately 7 a h/10 — 7 a fc/100. 

We can model the Doppler broadening by averaging 
over a Gaussian distribution of the one photon probe 
detuning A p and the two photon detuning, S = A p — A M . 
Written in terms of these two detunings (see appendix 
D for a full statement of the susceptibility including all 
nonzero detunings) the Doppler broadened susceptibility 
as a function of the mean probe detuning, A p , with mean 
pump detuning fixed at is given by 



1 



dA' p d6 x (1) (a;,5) 



-(A' p -A p f/(2*l) (S-A p f/(2v*) 



(43) 



group velocity index of n g 



- 14 400 in a Cs atomic va- 



Evaluating Eq. [43] numerically for a variety of vari- 
ances (see Figs. ITTJI and [TTj) . we can draw several conclu- 
sions about the constraints imposed by Doppler broad- 
ening. For one, the narrow features are unaffected by the 
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Two photon Doppler broadening 

to A" 




-0.10 



One photon Doppler broadening 

to x 



^^SSfZ^ 0?10 7ab 

M 




FIG. 10: (Color online) Here we have the two photon Doppler 
effect on the real and imaginary parts of the susceptibility 
near the narrow features. In both plots, we keep the one 
photon variance constant (to clarify the effect of the two pho- 
ton broadening) at <ta p = .00 17^6. The formatting of the 
lines is consistent between the two plots. In both, the blue 
dotted line corresponds to as = .00ly a b, the red dashed line 
has ag — .0057 a b, the brown dashed line has as = -0l7 a 6, 
and the solid green line (for which the features have van- 
ished altogether) has as = .057 a 6. For simplicity and gener- 
ality, we have used generalized populations rather than par- 
ticular models for open or closed pumping, with tys = -1 
and p c / c i — .8, and we have allowed the dephasings jc 
and 7c to vanish. The other parameters are fi M = 2^^, 

Q b = fi c = Jab/10. 



FIG. 11: (Color online) The one photon Doppler effect. 
Again, we keep the two photon broadening constant at as = 
.00l7a&. In this case we show a wide array of values for cta p , 
as the features are much less sensitive to the one photon 
broadening than the two photon broadening — indeed, each 
line represents variance an order of magnitude larger than 
the last. In both the top and bottom plots, the blue dot- 
ted line corresponds to <ta p = ■0l"fab, the red dashed line 
has a A p = -l^ab, the brown dashed line has <ta p = 7a&, and 
the solid green line has <ta p = 107 a !>. For simplicity and 
generality, we have once again used generalized populations 
rather than particular models for open or closed pumping, 
with ^Ps = .1 and p c i c i = .8, and we have allowed the de- 
phasings 7c and fc' to vanish. The other parameters are 
fi M = 2jab, &b — = Yab/10. 



single photon broadening, at least while <JA p /lab J$ 1- To 
see why this would be, consider the expansion of lm(x^) 
around ±fi&/2 in terms of the one and two photon de- 
tunings. We find, 

Im(x (1) (A p ,<S))«^(Re(^)-p^)x 

( 7a»OgA%+7C< 

^ { lab nyni + lc ,y + (5 T n b /2 + a*/n* (A p t n b /2))* 

(44) 

Under the assumptions that il^/fl^ <C 1 and A M — > 0, 
this expression reduces to the Lorentzian given in Eq. 1281 
In this case, Eq. HJdepends on <5 alone (because the term 
f2^/0^(Ap=pfif,/2) can be neglected in the denominator), 
which explains why the widths of the features are more 



sensitive to inhomogeneous broadening of the two pho- 
ton transition. Meanwhile, when a a p ^c / '^t approaches 
the width of the features, broadening of the single pho- 
ton process becomes important. Nevertheless, even for 
<7a /"fab ~ 10, our numerics show that the features will 
persist, albeit in suppressed form. 

The narrow features are far more sensitive to the two 
photon broadening. When as approaches the widths of 
the features, 7a6^c/^u+7c , they are rapidly broadened 
and ultimately washed out. One possible way of avoid- 
ing this difficulty is to work with a Doppler-free geometry, 
where the probe and the control lasers are copropagating. 
In this case, the two photon broadening vanishes and it 
would be possible to observe the (one photon broadened) 
narrow resonances even in a room temperature system, as 
per the discussion in the previous paragraph. It would 
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also be possible to work in a system with two photon 
Dopplcr broadening, so long as the temperature is suf- 
ficiently low. For fll/n^ w .01, it would be possible to 
work in a laser cooled system; for smaller f2 c /f2p, a quan- 
tum degenerate system would be necessary to realize the 
gain. 

The anomalous dispersion, meanwhile, is more robust 
under Doppler broadening. Just as when one increases 
the homogeneous broadening of the narrow lines, the 
slope of the real part of the susceptibility in the anoma- 
lous regime is unchanged even for large inhomogeneous 
broadening. It is only when the features are washed 
out entirely by the inhomogeneous broadening that the 
anomalous dispersion is lost. In particular, the anoma- 
lous dispersion is not reduced when the variance in the 
one photon detuning, ua p , is as large as 107 a fc. The dis- 
persion is more sensitive to the two photon detuning, but 
only because the features themselves vanish entirely for 
a smaller value of as- This means that the anomalous 
dispersive regime should be readily observable and un- 
suppressed even in a room temperature gas, provided one 
uses a Doppler-free geometry to eliminate the two pho- 
ton broadening. Our numerical calculation predicts that 
even in the absence of co-propagating lasers, the disper- 
sion will be unchanged for as < 7 a b/100, even though the 
resonances will be severely broadened. Thus it should be 
possible to observe the anomalous dispersion in a laser- 
cooled gas without a Doppler-free geometry. 



VI. SUMMARY AND CONCLUSIONS 

As a final note on experimental realizability, we pro- 
pose several possible level structures. The levels we pro- 
pose here are for 87 Rb, though the equivalent levels in 
23 Na would work equally well. The first proposal is a 
modification of the structure used in the EIT experiment 
in Ref. (43[, which we used above to make numerical 
predictions for the values of the negative group delay in 
our system. For the closed pumping system, one might 
use Di line with \a) = \5P 1/2 ,F = 2,m f = 2), \b) = 
\5S 1/2 ,F = l,m f = 1), \V) = |55 1/2 ,F = l,m f = 0), 
|c) = |55 1/2 ,F = 2,77^ = 1), and |c') = |5S 1/2 ,F = 
2,m/ = 2). Alternatively, one might use the D 2 line 
with \a) = \5P 3 /2,F — 2, m/ = 2) and the same ground 
states from the 5Si/ 2 manifold as just given for the D± 
line. Then \a) decays only to \b), |c), and |c'), as de- 
scribed by our model. (Note that choosing levels where 
| a) decays also to \b') will not change qualitatively any 
of our results provided fib ^ and A;, = 0.) For both 
of the configurations we propose, both the probe and the 
control beam would have to be <r+ polarized. 

In conclusion, we have studied the effects of various 
pumping configurations on the linear response of a driven 
five-level atom. We have found that when the population 
of one of the ground states, |c'), becomes large, it pro- 
duces two amplification resonances without population 
inversion in the bare state basis. We analyzed the depen- 



dence of the population of |c') on two different pumping 
configurations. Moreover, we have shown that in the re- 
gion between the two gain lines, it is possible to tune the 
system to permit anomalous dispersion. We have stud- 
ied the effects of Doppler broadening on this system, and 
concluded that the anomalous dispersion is robust under 
both homogeneous and inhomogeneous broadening of the 
gain lines, so long as the gain lines do not vanish. 

C.P.S. acknowledges support for this work from the 
National Science Foundation award no. 0757933. 



APPENDIX A 



The equations of motion for the system described in 
section [TT] are given by, 



.dpaa 

.dpbb 
. dp b 'b> 
.dpec 



^{Pca ~ Pac) + ^-(Pab ~ Pba) 



~^{PbV ~ Pb'b) ~ 

— i^\Pbb' - Pb'b) 



(Pab - Pba) 



a 



dt 

.dpc>c> 



{Pc 



a. 



Pc 



{Pc 



a 



{fide - Pec') 



(Ala) 
(Alb) 
(Ale) 
(Aid) 
(Ale) 



and then the off-diagonals, 



.Opab 
. dpab' 

''~m~ 

.dpea 
. dpc'a 

'~dt~ 
.dpcb 



ApPab Y Pcb + ~P ab ' + — {Paa - Pbb) 



~2 



2 



(A p - A b )pab' TTPcb' 



a, 



\ - '*V ~ _ ^jp l ~ _ ~ \ _ ^£ - 

~t-^fj,Pca i 2 Pcb 2 \P aa Pec) 2 Pc'a 



(A c - A^J/Va + ~ Pc ' c + ~2 Pc ' b 2~ Pca 



(Alf) 
(Alg) 

Pc'a 

(Alh) 
(Ali) 



/ A A \ ~ . ^b ~ ~ QfJ, ~ 

(Ap - Af,)p cb + —Pcb' + —Pea ~ -^Pab 

Mc s 
- ~^Pc'b (Alj) 
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dp, 



dt 



= (A p - Afe - A^pcb* + -^Pcb — ^-Pab> 



.dp, 



c'b 



dt 



.dp 



'c'b' 



dt 



.dp b b' 

l -dT 

. dpc'c 



-jPc'b' 



(Alk) 



(A p + A c - A^pc'b + ^-Pc'b' + ^Pc'a 



c ~ 
-Pcb 



= (Ap - A b + A c - A^)p db , + -^-Pcb 



Y Pcb ' 



(All) 



(Aim) 



= -A b p w - ^pab' + ^{Pbb - Pb'b') (Aln) 

= A c p c > c + ^-Pc'a - ^-(Pcc - Pc'c') (Alo) 



subspaces of the Hamiltonian via the matrix 



D 



where 

















cos 6 b 


sin 6 b 








— sin 6 b 


cos 6b 














cos 6 


\o 








— sin 



\ 



sin# 



(Bl) 



cos^i±^ sin^-/ 1 -^ 



off 



We have defined detunings A p = to a — u> b — v p , A M = 
Lu a - Lo c - Vp, Ab = uJb> — uj b — v b , and A c = u c < —u> c — v c . 



APPENDIX B 



and 



n 



f = V / A|7nf. 



We transform the equations of motion to first order 
in fip/Q^ by diagonalizing the {|6),|6')} and {|c),|c')} 



In this basis, the Hamiltonian becomes 



DHD^ = 



-Q p cos 6b 



( 2cu a 

Q p cos 6b 2uJb + Ab + 2v p — il 



y fi M sin0 c 



fip sin 6b 
-fi M cos 6 r 



cir 

b 



dp sin 6b 




—Clfi cos( 





2uj h + A b + 2v v + Qf 
2uj c + A C + 





2^ - nf 



£l M sin6> c \ 






2uj c + A C + 2u fJt + Of J 
(B2) 

Written, for now, without decay (which will require some approximations to transfer into this basis), we find 

= *- Ap ~~ ~T + ~2~^ aB ~ -f( cose bf>BB - smdbPB'B) - ^(cos6» cj 5 CB - sm6 c p c > B ) + y cos 6> fc/ 5 aa (B3a) 



.dpcB 



= (-A„ + Ap 



A b Qf a c Qf n p n„ 

— + 2 ~2 2~' PCB 2 cos0bPCa ~ 2 cos9cPaB 



.dp^B = + Ap _Ab + ^ + Ac + ^_ )p>s + ^ cos6l 6 pc'« + ^ sin0 c /5 aB 



at 



2 



Sip 

~2~ 



Op 
2 



(B3b) 
(B3c) 



and 



at v p 2 



■ dpcB' / a i a 

= (-A„ + A p 



Qeff Q ^ Q 

-^-)paB - -^(cos6» b( 5 B B' - sin6» b p B / B /) ^-(cos6 c p CB ' -sm9 c p~c'B>) - -y sin6> b/ 5 a 



ft? ff A, OS ff 



+ 



SI 7) 



a. 



,0Pcb' t a i a A 6 fif* A. 



at 



(-A„ + Ap 



->cff 

2 



2 



£7p 



— - - ■ ^ sin 6> fc/ 5 C / a + ^ sin 6» c p aB - 



(B4a) 
(B4b) 

(B4c) 



In this partially diagonalized basis, we have p a b — cos 6bp a B — sin 6bp aB ' , pea = cos 6 c p ca + sin 6 c p c i ai and pc'a = 
cos c /? c 'a — sin 6 c p ca . Meanwhile, 

( Pbb Pbb'\ _ /cos 2 6 b pbb + cos 6 b sin 6 b (pb'b + Pw) + sin 2 6 h pb'v svn 6 h cos 6 b{pb'v - Pbb) + cos 2 6 h pw - sin 2 6 b pb'b 
\Pb'b Pb'B'J \sind h cosd b (pb'b' - Pbb) + cos 2 O b Pb'b ~ sin 2 Q b p w cos 2 6 h Pb'V ~ cosd b sm6 b (p b > b + p w ) + sin 2 6 b Pbbj 
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APPENDIX C 



APPENDIX D 



To incorporate decay in the dressed basis, we have 
to make certain assumptions, as described in the text. 
These amount to saying that j ab « J a b'- and that 



7h ~ 76') lie ~ 7b' C : and 7b c ' ~ Hb'c' so that we can 



take 7 c6 w -f cb , = -f C and 7 c / fc w 7^/ = 7c- Under 
these approximations, the contributions from decay and 
dephasing are 



For completeness, we present here a complete and gen- 
eral solution for the reduced susceptibility, , including 
arbitrary detunings of all fields. We find 



ipaB ' 


^ -ilabPaB 


(Cla) 


ipaB' ' 


^ -ilabPaB' 


(Clb) 


ipCB ' 


- — i(7o cos 2 9 C + 7c sin 2 9 c )p C B 






- i(lc> - 7c) cos 6 C sin 9 c p c >B 


(Clc) 


ipc'B ' 


v> -i(j c sin 2 C + 7c" cos 2 9 c )p c >B 






- i (7c - 7c ) cos C sin c pcs 


(Cld) 


ipcB' ' 


- -i(7c cos 2 C + 7c sin 2 9 c )p C B> 






- i (lc> - 7c ) cos C sin 9 c p c >B' 


(Cle) 


ipCB' ' 


- — i(7c sin 2 C + 7c cos 2 9 c )p c >B> 






- i(lc> - 7c) cos 9 C sin 9 c pcB' 


(Clf) 



2"fab 
{ In 



(cost 



% b - sin f 



(Dl) 



Here. 



and 



PaB = ((sin 9 b p B >B - cos 9 b p BB ) (- (2i lc + A b - A c - 26 - Qf) (2ij c > +A b -A c -26- Qf) 

+2i cos(20 c ) ( 7c - 7 c) Qf + (Qf) 2 ) + cos 9 b cos 9 c p Ca (2«7c< + A b - A c - 25 - Qf - Qf) n, 
- cos0 b sin 9 c p C 'a (2*70' + A b - A c - 2«S - nf + nf) fi p ) 



n 



Y~ ((sm9 b p B >B> 



25 + nf ) (2i 7 c + A fc 



5 + n? 



2^ + n? 



nf) n, 



os6 b pBB') (- (2i 7c + A 6 - A 

+2z cos(20 c ) ( 7c - 7 c) nf + (nf ) 2 ) - sin b cos 9 c ~p Ca (2i lc , + A b 
+ sin 9 b sin c p o , o (2*70* + A b - A c - 26 + Qf + Qf) 

where now 

Z± = sin 2 9 C (2i 7 c + A b - A c - 2<J T nf + Qf) n 2 + 2z cos 2 9 C sin 2 C ( 7c - 7c) (2t (7c - 7c) 

x (2 7ab + A b - 2A P t nf) - n 2 ) + (-2i sin 2 c7c - 2 Z cos 2 c7 c - A b + A c + 2(5 ± nf + Qf) 
x ((2i 7ah + A b - 2A P t nf ) (2i cos 2 9 clc + 2i sin 2 c7 c + A fa - A c - 26 t Qf + Qf) ~ cos 2 9 C Q 
We have defined the two photon detuning, 6 — A p — A M = lu c — uj b + v^ — v p . 



(D2) 



(D3) 



(D4) 



In the special case that cos9 b psB — svn9 b pB'B — sm9 b pB'B' + cos9 b pBB' , we can define tys = \/2(cos9 b pBB — 



sm9 b p~B' b)- Likewise, when sin 9 c pc a = cos0 c pc a , we can define tyc 
on the relatively simple form familiar from the body of the paper, 

n, 



2S1 



sin 9 c p Ca ■ Then Eqs. |D2] and [D3] take 



Pas = (Vb ((2i7c + A b - A c - 26 - Qf) (2i lc - +A b -A c -25- Qf) - 2i cos(20 c ) ( 7c - 7c) Qf) 



(Qf ) 2 (cos 9 b %c- < VB)) 



(D5) 



and 



o aB > = (!Pb ((2 t7 c + A b - A c - 2<5 + nf) ^ 7c , + A b - A c - 6 + Qf) - 2i cos(20 c ) ( lc - yc>) Qf) 

■\efi\2 



+ (Qf) 2 (sm9 b ^ c -VB)) 



(D6) 
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